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Abstract 

Variational calculations using Gaussian wave functionals combined with an ap- 
proximate projection on gauge invariant states are presented. We find that the 
energy exhibits a minimum for a wave functional centered around a non vanish- 
ing background magnetic field. We show that divergences can be removed by a 
renormalization of the coupling constant. The resulting expectation value of the 
gluon condensate is found to be in qualitative agreement with phenomenological 
estimates. 
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The application of the variational method to gauge theories is generally plagued 
by the difficulty to implement in a calculable way the requirement of gauge invariance of 
physical states [Q, |T], In the present letter we propose to use an approximate projec- 
tion method developped in 1962 by Thouless and Valatin 0. This method deals with 
the restoration of rotational invariance when deformed solutions are obtained in nuclear 
Hartree-Fock calculations. 

In the case of gauge theories this formalism allows one to explore a class of invari- 
ant trial vacuum functionals for which closed expressions can be obtained. This class is 
generated by projecting Gaussian wave functionals onto the subspace of gauge invariant 
states. This class of trial functionals has already been explored by Kogan and Kovner 0], 
and Brown and Kogan Q . Since these authors use different approximation schemes their 
approach and ours appear complementary. 

We consider the functional Schrodinger description of the SU(A^) Yang-Mills theory. 
In the temporal gauge ^4° = 0, the canonical coordinates are the vector potentials v4^(x) 
and the electric fields i?*(x), which we shall often write as color matrices in the adjoint 
representation: e.g., A* = A^T'^ (the color indices a, b, . . . run from 1 to M = A^^ — 1). 
The generators T" of the color group are taken to be Hermitian and traceless; they satisfy 

[T", T''] = ir^'^T", Tr (T"T^) = W^ (T")b, = -1^'''. (1) 

The Hamiltonian density reads 

n{^) = ^{E^Eli^) + BlBli^)}, (2) 

with the color magnetic fields = and F'J = d'Ai - d^A^ + gfabcK^c 

{g denotes the coupling constant). By also introducing the covariant derivative = 
9* — igA^, we can rewrite the previous formulae in matrix form: e.g., F'-' = F^^T"- = 
[D\D^]/{'ig). 

In the Schrodinger representation, the states are represented by functionals of Aj,(x), 
^[A], and the electric field is acting on such states by functional differentiation: 

=^^^^[A]. (3) 

The hamiltonian H commutes with the generator Q of time-independent gauge transfor- 
mations, 

g{^) = V-E(x) + zg[A,E% (4) 

so it is possible to diagonalize H and Q simultaneously. We are here interested only in 
those eigenstates \I' of which are gauge invariant, i.e., 

^(x)vl>[A] = 0, (5) 
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or in superpositions of them, to be referred to as physical states. 
The variational principle consists in the inequality 

which holds for any functional ^[A] from the physical Hilbert space. In this equation, 
E^ac is the ground state energy, assumed to be non-degenerate. In practice, however, one 
has to restrict oneself to Gaussian wave functionals, the only ones which allow for the 
analytical computation of (H). These have the form 

v^[A] = exp |- i I d'xd'y [A^(x) - A^(x)] {G-')"^' (x, y) [Ai{y) - Ai{y)] } , (7) 

where the background field Al^('x) and the kernel G^^ (with matrix elements (G^^)"* (x, y)) 
are the variational parameters. 

The expectation value of the Hamiltonian density in the Gaussian state is 
{^\H{^)\^) = |B- B(x) + |Tr(x|G'-i|x) + iTr(x|irG'|x) 

+ sl (Tr{5,T'^(x|G|x)})^ (8) 

+ ^Tr {5*T"(x|G|x)5^T'*(x|G|x)}. 

In this equation B is the magnetic field associated to the center A and is the spin one 
matrix whose elements (j, k) are given by iSijk- The notation Tr in equation (j^) implies 
a summation over both the color and the spacial indices. For instance, Tr(x|G"^|x) = 
Si,a (G~^)ff (x, x). Finally, the operator K is the second derivative of the classical energy 
with respect to the center A"}. It is given by 

K = {-iS ■ X>)2 - S ■ B, (9) 

where 

V' = d' - igA' (10) 

denotes the covariant derivative defined by the background field A^ = A^^T"'. 

In the case of non-Abelian gauge theories, however, the Gaussian functionals suffer 
from a major drawback: they do not satisfy the requirement of gauge invariance (^). It 
is thus necessary to project them onto the subspace of gauge invariant states. This is 
achieved by means of the formula 



^o[A] = -^/l?[f/(x)]vI/^[A], (11) 
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where the functional integration is performed over the unitary N x N matrix field f/(x), 
with the adequate group invariant measure, and A/" is a normalization factor. The inte- 
grand in eq. (|rT|) is the gauge-transform of ^[A]: 

^c/[A] = ^[UAU+ + -UV U+]. 

The expectation value Ep of the energy in the projected state is given by the following 
formula 

_ sv[u{^)]<m\^u> ...^ 

which should replace eq. (H) in practical calculations. Unfortunately, eq. (|l^) cannot be 
evaluated in closed form because the functional integral over the group is not Gaussian 

i- 



In what follows, we shall propose an approximation to eq. ([I2[) which is inspired 
from techniques used in nuclear physics to calculate the zero point rotational energy of 
deformed nuclei. The starting point is the formula for the collective rotational energy of 
a classical rotating rigid body 

E = —lijUJiUij, 

where Tij is the inertia tensor and Ui the angular velocity along the i-th axis. When de- 
formed solutions are obtained in solving the nuclear Hartree-Fock equations, the Thouless- 
Valatin formalism provides a quantum extension of this formula. It gives the energy gain 
when projecting on zero angular momentum states as |^ ^ 

, ^ < > < Jy > < J"? > 

AE = — I — + — ^ — + — 4 — . (13) 

In this formula is the angular momentum operator along the x axis and 2^ the corre- 
sponding moment of inertia. This quantity is the polarizability 

X. = lim^^ (14) 



of the system in the presence of an external constraint uJxJx- Equation (13) is an ap- 
proximate expression valid in the limit of a large fluctuation in the angular momentum 
< J2 » 1 §. 

In the case of gauge theories it is straightforward to generalize the above formula 
(see for more details). The moment of inertia is now a matrix in color and position in 
space defined as the polarization tensor 

r'U y) - ^^^"(^^^ 1 fl5) 
^ Su^iy) '-=0' ^^^^ 
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in the presence of an external constraint 

H,.t = J d'y AyW\y)- (16) 

In the equations above, Q is the generator of time -independent gauge transformations (cf. 
eq. (^). The gain in energy when projecting a wave functional ^E'[A] onto the subspace 
of gauge invariant states reads, in analogy to eq. (]T3|), 

AE = J d'xd'y{^\gy^)g\ym{a, x|^|6, y). (17) 

The expectation values in eqs. (|15]) and (|l^ are all Gaussian and can be exactly computed. 
Nevertheless, the ensuing calculation of the moment of inertia is extremely tedious, which 
invites us to look for a convenient approximation. Namely, we shall assume in what 
follows that the condensate is large enough so that we can replace the field operator A by 
its mean value A in the expression (|^) for the gauge generator. We thus define: 

^(x) = V-E(x) + ig[A\E']. (18) 

It is then straightforward [0] to obtain the corresponding, approximate form of the inertial 
moment: 

(a, x|J|6, y) ^ ^^^5(y)^^ Uo = («' ^|n>, y) (19) 

where 11^ denotes the square of the kinetic momentum Ilj = iT>j = idj+gAj. Furthermore, 
from eq. (^ and (0), one readily obtains: 

{^m^)g''{y)\^) = i (a, x|n,Gr.in,|6, y) , (20) 



so that the projection energy (^) can be written as: 

1 



AE = - I d^x (a, x| ( Ui—Uj GT^^\a, x) . (21) 



The Thouless-Valatin correction thus changes the term Tr(l/8(j) of equation (P) into 
Tr(g/8G') where 

Qij = Sij - Ui—Uj. (22) 

It is easily checked that Q^=Q and that the action of Q on a transverse vector field v4"(x) 
(i.e., such that IliAi=0) gives back A'}. We will refer to Q as projector on transverse 
fields. We conclude that the Thouless-Valatin correction makes the kinetical part of the 
energy transverse with respect to the mean-field covariant derivative Dj. 

Further note that the operator K is transverse as well. In fact, eq. (^ can be 
rewritten as 

K,, = U%, - UiU, + 2[Ui,U,] , (23) 
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which is manifestly transverse: HiKij = 0. Then, by minimizing the projected energy 
E — AE with respect to G, we immediately find that the kernel must be transverse, 

V'GT.^ = 0. (24) 

This condition has a simple interpretation: it is equivalent to ^(x)^[A] = 0, showing 
that the Gauss law (||) is approximately satisfied by the Gaussian (|^). 

Eq. (1^) requires the operator G^^ to be a functional of the background field Ai. A 
specific Ansatz for G~^ can be obtained from the variational equation for the projected 
energy E — AE. To this aim, we write G^^ = QG^^Q (and similarly Gt = QGQ), 
and consider the gap equation for Gt- In order to exploit this equation, we consider a 
restricted variational space defined by the following background field: 

A^ = 0, Ay = x T^B, A, = 0. (25) 

This corresponds to a constant magnetic field in the z-direction and in the third color 
with a strength B which is a variational parameter. For such a background, one can show 
1^ that the most general kernel Gt which is consistent with the gap equation is of the 
form 

= K + PQS - BQ = Q (U^ + aS ■ B)Q, (26) 

where the dimensionless coefficient (3 = 2+a is a new variational parameter. The accept- 
able range of values of /3 is > 1 . Indeed, for /3 < 1 the previous equation is undefined 
because of the occurrence of negative modes in the right hand side P, p!0 . 



The expression of the energy density involves quantities such as < xIGtIx > and 
< xlGj^^lx >. A proper definition of these quantities requires a regularization scheme. 
A convenient one in the present context is Schwinger's proper time representation which 



gives 



(xIGtIx) = ^ r ^(x|e-*(™«-^«)|x), (27) 



and 



IT Jl/A2 



Integrations for small values of the proper time t have been regularized by the introduction 
of a cutoff A. In the limit of a large cutoff one derives 

Tr(G?i) = ^ |2A^ - c;v(y -2(3 + ^)g'B' In ^ + • • j , (29) 
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and 

(x|G^|x) = I^A' + (1 - f )S ■ B In ^ + (x|G^|x) + . (30) 

In this last equation Gi? is the finite part of Gt in the hmit of a large cutoff. An ap- 
proximate closed formula for this term can be obtained by ignoring the projection on 
transverse states in the Ansatz (^6]) for the kernel of the Gaussian. In this case one finds 
indeed 0] 

The resulting energy density reads 



+ 7dw4c%B'{l - f)Vln' f + (^) WA2Tr(x|G^|x), 



(32) 



(647r2)2'-2-Af-^ V-^ 3 y y B ' V 32,72 

where Cn= Tr(T(3)r(3)) = N for the group SU(N). Ignoring the -B-independent divergent 
piece ~ A^, the remaining expression in eq. (^) should become finite after the renormal- 
ization of the coupling constant: that is, (\E'|7i(x)|\I') must become finite after sending the 
bare coupling strength to zero logarithmically in the limit of a large cutoff. But obviously 
this procedure cannot work for the last term in eq. (|3^) , which involves a quadratic ul- 
traviolet divergence. We conclude that the renormalizability of (\E'|?i(x)|\E') requires the 
trace of the finite part of the kernel to vanish identically, which then gives a condition 
on the variational parameter f3= 2 + a. Using the approximate expression (0), this 
condition reads: 

rds( eBs ^2 cosh aeBs)- 3} = 0. (33) 
JO I sinh eBs J 

An approximate solution of this equation (with non-vanishing B) can be found by ex- 
panding the hyperbolic functions to second order. One thus obtains: 

a^--, i.e. /3~2±^. (34) 
2' ^ V2 

Remarkably, these values of (3 guarantee that, in the present calculation, there is no 
problem with the negative modes (recall that, for (3 = 0, one faces infrared divergences 



at large proper times [|T3[; see, e.g., eq. (plD ). The minimum in the energy density can be 
seen to correspond to the value /? = 2 — 1/ v^. For this value of /?, the energy density is, 
in the case of the SU(3) gauge group. 
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where only the -B-dependent terms have been kept. The minimum in the energy density 
occurs for B = Bmin with 



Bmin = exp( Xmin), (36) 

9^C3 



and 



__9 

The value of the energy at the minimum is 
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X™„ ~ 0.363 + + Oig^). (37) 



< n >rmn= ^ {(^ + V^)X™„ " f " V^} < 0- (38) 

In order to obtain a finite value of the magnetic field at the minimum when the cutoff 
A is sent to infinity, namely Bmin= / -\/e, the coupling constant g must run with the 
cutoff according to 

^ ^ In— = 0.363. (39) 



167r2 /i 

We can now define the effective coupling constant for a given scale /i by looking at the 
coefficient of B^ in the expression of the energy density and by writing the quantity 
ln(AV5) as ln(AV/i2) + ln(/iV5). The result is 

^ - ^ ^ (13 + 4v^)c3ln^ + ^^^|^(3 + 2V2)c^ln2^. (40) 



From this equation we obtain the beta function 

%R 3 1 



l3{gM) = /i^ = -gl-^csV52m, (41) 



in the limit g\/ g"^ 0. It differs from the perturbative one loop result 

IMS) = (42) 

This is not unexpected in view of the fact that our vacuum state exhibits a non vanishing 
value of the condensate and a kernel which largely differs from the perturbative one. In 
fact the perturbative result can be recovered in our approach by simply ignoring the term 
PSzT^^^ in the trial Gaussian kernel. This is however not legitimate in the variational 
context because (3=0 belongs to a region where the kernel exhibits unstable modes. Our 
result also differs from that of Kogan and Kovner who use a variational Ansatz which 
does not include the possibility of a condensate. 
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By looking at the energy expression at the minimum we find that the quantity 
< Ti >min / 9^ has a hmit 

^ ^ f '"^ ~ -4.2 X lO'V (43) 
9^ 

when the cutoff goes to infinity, i.e. is a renormahzation group invariant quantity. We 
can compare this result to the value of the energy density in the bag model e = —B 
where B is now the bag constant With 5= (240 MeV)^, the bag model leads to e= 



-3.3 X 10-3 GeV^ This gives /i=0.941 GeV. 

iFiom our approximate solution for the vacuum functional, it is easy to calculate 
the value of the gluon condensate: 

< F;,F;^ >= -2 < E^E^ - BfB^ > . (44) 

At the minimum, we obtain for this quantity a finite value 

< F^^F;^ >= + V2)5Ln^™n, (45) 

which gives 

< >^ V, (46) 

where /i ~ 1 GeV. This result is in qualitative agreement with the phenomenological value 
0.5 GeV^ of the gluon condensate . 
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